Magnetization plateau in the spin ladder with the four-spin exchange 
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The magnetization process of the 5=1/2 antiferromagnetic spin ladder with the four-spin cyclic 
exchange interaction at T = is studied by the exact diagonalization of finite clusters and size 
scaling analyses. It is found that a magnetization plateau appears at half the saturation value if 
the ratio of the four- and two-spin exchange coupling constants J4 is larger than the critical value 
Jic =0.05±0.04. The phase transition with respect to J4 at Ji c is revealed to be the Kosterlitz- 
Thouless-type. 
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The field-induced spin gap is one of recent interesting 
topics on the one-dimensional (ID) quantum spin sys- 
tems. The gap can be detected as a plateau of the mag- 
netization curve at low temperatures. The appearance of 
such a magnetization plateau was theoretically predicted 
in several systems; the anisotropic S = | antiferromag- 
netic chain , the S — 1 bond-alternating chain 0] , the 
S = h ferromagnetic-ferromagnetic-antiferromagnetic 
chain Q, the frustrated bond-alternating chain H, the 
three-leg ladder Q , and the frustrated two-leg ladder . 
In particular the two-leg ladder attracts a great interest 
in the context of the superconductivity in a doped system 
H . The standard S — \ uniform antiferromagnetic spin 
ladder has the spin gap of the lowest excitation from the 
nonmagnetic ground state (GS), which leads to a plateau 
at the magnetization m = p-pd|. A strong coupling 
approach Q indicated an additional plateau at half the 
saturated magnetization due to the next-nearest antifer- 
romagnetic coupling which yields the frustration. In this 
paper we show another possibility of the magnetization 
plateau in the two-leg spin ladder, which is induced by a 
four-spin exchange interaction. 

A four-spin exchange interaction described by a prod- 
uct of two-spin exchanges in a spin ladder was investi- 
gated by a field theoretical approach fl^| . It indicated 
the possibility of a different type of massive phase from 
the Haldane phase jl3| in the nonmagnetic GS, but the 
state in a strong magnetic field was not discussed. On 
the other hand a mean field analysis [|l4| suggested that 
the S = h triangular lattice antiferromagnet would have 
a magnetization plateau at half the saturated magneti- 
zation, if there exists a four-spin cyclic exchange inter- 
action. It was verified by the exact diagonalization [fl5|| . 
The recent experiments revealed that such multiple-spin 
exchange interactions are realized in the two-dimensional 

(2D) solid 3 He @,0 and thc 2D w W eT solid of clcc " 
trons formed in a Si inversion layer JL8[ , as well as the 
bec 3 He [[l9|. In order to test the possibility of a similar 
magnetization plateau in ID quantum spin systems, we 
consider the S = i uniform antiferromagnetic spin lad- 
der with the four-spin cyclic exchange at every plaquette. 
The magnetization process of the system is described by 



the Hamiltonian 
H = Ho + Hz , 

Hq = V^(Sij ■ Sij+i + ■ S2J+1 + Sij ■ S2,j) 
3 

■■i.yji\., ■ (i) 

H z = -ifyi(^i,j + Slj), 

3 

where P^j is the cyclic permutation operator which ex- 
changes the four spins around the j-th plaquette as 
Si j — > Sij+i — > S2J+1 — > — ► Si j, J4 is the 
strength of the four-spin exchange and H is the applied 
magnetic field normalized by the two-spin exchange cou- 
pling constant. We assume J4 is positive, as it is in the 
solid 3 He. This system subjected to the periodic bound- 
ary condition is studied by the exact diagonalization of 
the finite clusters and the size scaling of the low-lying 
energy spectra. For L x 2-spin systems, the lowest en- 
ergy of Ho in the subspace where . (Sfj + Sfj) = M is 
denoted as E(L,M). Using Lanczos' algorithm, we cal- 
culated E(L, M) (M = 0, 1, 2, • • • , L) for even-L systems 
up to L — 16. The macroscopic magnetization is defined 
as to = 

The nonmagnetic GS of the system ([!]) with J4 = 
is in a massive phase equivalent to the Haldane phase of 
the S = 1 antiferromagnetic chain and the low-lying ex- 
citation has a finite energy gap for to = 0. On the other 
hand, the magnetic GS is always gapless [^0|,^l) except 
for the saturation. Thus the magnetization curve has a 
plateau at m = 0, while no other plateau appear up to 
to = 1, as far as J4 = 0. The four-spin exchange, how- 
ever, is expected to induce a plateau at to = s, because 
the interaction stabilizes the state, mentioned in 

Ref. |TJ] , of the four spins around every plaquette within 
a mean field argument. We concentrate on the plateau 
at m = i, rather than the nonmagnetic GS. 

The magnetic excitation gap giving 8M = ±1 of the 
L x 2-spin systems described by the total Hamiltonian H 
is given by 



1 



A± = E(L, Mil)- E(L, M) t H. 



(2) 



The sum A 



A_ is a good order parame- 



For the gapless system in the thermodynamic limit, the 
conformal field theory |22| (CFT) predicts the asymptotic 
form of the size dependence of the gap as A± ~ 0(1/ L) 
with fixed m — M/L. When H + and H- are defined as 



E(L, M + 1) - E(L, M) -> H+ (L 
E(L, M) — E(L, M — 1) — > H- (L 



oo), 



(3) 



H+ and LT_ has the same value and it gives the magnetic 
field H for the magnetization m in the thermodynamic 
limit. In contrast to the gapless case, if the system has a 
finite gap even in the infinite length limit, A+ and A_ are 
still finite for L — > oo. It leads to the difference between 
H+ and H- and a plateau appears for H- < H < H + 
at m = M/L in the magnetization curve at T = 0. 
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FIG. 1. Scaled gap LA versus the strength of the four-spin 
exchange interaction J4. 
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FIG. 2. L-dependent fixed point Ja c l,l+2 of the gap for 
SM = ±1 (circles) and SM — ±2 (squares) are plotted versus 
1/(1/ + 1) to determine J^ c in the thermodynamic limit. The 
estimated value is J^ c — 0.05 ±0.01 for SM = ±1, which does 
not well agree with the result for SM = ±2 J 4c = 0.01 ± 0.01. 
Thus we conclude J 4c = 0.05 ± 0.04. 



ter to investigate the plateau-nonplateau transition with 
the finite-size scaling M, because A corresponds to the 
length of the plateau in the magnetization curve in the 
thermodynamic limit. The scaled gap LA of finite sys- 
tems = 8 ~ 16) at m = 1/2 is plotted versus J4 in 
Fig. For J4 > 0.2 the scaled gap obviously increases 
with increasing L, which means that a finite gap exists 
in the thermodynamic limit. For small J4 around the 
region < J4 < 0.1, the scaled gap looks almost inde- 
pendent of L. It implies that the system is gapless at 
a finite region of the parameter J4, which is reminiscent 
of the Kosterlitz-Thouless (KT) transition [M . Accord- 
ing to our precise analysis, the LA curves for L, and 
L ± 2 have an intersection in the region < J4 < 0.1 
for each L. Thus the critical point J 4c can be estimated 
by the phenomenological renormalization group (PRG) 
equation 24 



(L + 2)A L+2 (J 4 ) = LA L (J 4 ). 



(4) 



We define JacL.l+2 as the L-dependent fixed point of 
(U) and it is extrapolated to the thermodynamic limit. 
JacL,l+2 is plotted versus 1/(L + 1) as solid circles in 
Fig. ||. Although the convergence of Ja c l,l+2 with in- 
creasing L is not good, the least square fitting of the form 
JicL,L+2 ~ Jic + A/(L + 1) gives the extrapolated result 
J4 C = 0.05 ±0.01 as the dashed line in Fig. g. To test the 
precision of the value, we did the same analysis using the 
gap for SM — ±2 instead of A± as solid squares shown 
in Fig. ^| where the fixed point can be obtained only 
for L > 10. It gave J 4c = 0.01 ± 0.01 which is not well 
coincide with the above result, which implies that the 
available system size is not enough to determine J4 C with 
the fitting of 1/(L ± 1). (Such a difficulty of the precise 
decision of the critical point by PRG is sometimes due 
to the logarithmic size correction in the case of the KT 
transition. pSj ) Assuming that the system is gapless for 
J4 = 0, we conclude J4 C = 0.05 ± 0.04 within the present 
analysis. 

We present the GS magnetization curve in the ther- 
modynamic limit for J 4 =0 and 0.1. In the latter case 
the magnetization plateau should appear at m — | in 
contrast to the former, as discussed above. Note that 
the four-spin exchange interaction reduces the spin gap 
just above the nonmagnetic GS. According to our present 
analysis based on PRG, the gap for m — vanishes at 
a critical value J4 C , which should be distinguished from 
J4 C for m = i, and the nonmagnetic GS will belong 
to a different massive phase from the Haldane phase for 
J4 > J4 C . The critical value J4 C , however, is obviously 
larger than 0.1. Thus even in the case of J4 = 0.1 the 
spin gap due to the Haldane mechanism still exists for 
m = and we can use the same method to give the mag- 
netization curve as used for the S = 1 antifcrromagnetic 
chain §|. 

For J4 = 0.1 the left hand sides of the form (||) calcu- 



lated for m = 



I and 



are plotted versus 1/L 



2 



in Fig. |[ It shows H + = H_ except for m = and i. 
Thus we take the mean value of the two for the magnetic 
field H for each m. 
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FIG. 3. E{L,M+1)-E(L,M) and E{L,M)-E(L,M-1) 
plotted versus 1/L with fixed m for J 4 = 0.1. The dashed 
curves are guides to the eye. The extrapolated points for 
m — 0, m = 1/2— and m = 1/2+ correspond to the 
results of the Shanks' transformation H c i = 0.15 ± 0.03, 
H- = 1.84 ± 0.06 and H+ = 1.99 ± 0.09, respectively. 

Since the nonmagnetic GS is massive for J4 =0 and 0.1, 
the size correction of H+ decays faster than j- as shown 
111 Fig. § Thus we use the Shanks' transformation 1291 
P' n = (P n _iP n+ i - P„ 2 )/(iVi + P n +i - 2P n ) twice for 
the sequence E{L, 1) - E(L,0) for L =6,8,10,12 and 14, 
and obtain H cl = 0.503 ± 0.003 and 0.15 ± 0.03 for J 4 =0 
and 0.1, respectively. The saturation field H C 2 is given 
by the L-independent quantity E{L, L) — E(L, L — 1). 
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FIG. 4. GS magnetization curves in the thermodynamic 
limit for J4 = and 0.1. The latter has the magnetization 
plateau at half the saturation value. 

In the case of J4 = 0.1, for m = 1/2 H + and £/_ are 



obviously different and the size correction decays faster 
than 1/L, as shown in Fig. [j], which is consistent with a 
finite gap. Then we estimate H + and H_ by the Shanks' 
transformation and get H + = 1.99 ± 0.09 and if_ = 
1.84 ± 0.06. For J 4 = H + and i?_ correspond even at 
m = 1/2. We present the results for J 4 =0 and 0.1 in 
Fig. [|, where we also used the values of H for m =1/6, 
3/8, 5/8, 5/6 and 7/8 which are estimated by the same 
method as mentioned above. The curve has a plateau at 
m = 1/2 (H-<H< H+) for J 4 =0.1. 

Our present PRG analysis shows that the the gap does 
not behaves as A ~ (J4 — J 4c ) iy . If we define the the 
size-dependent exponent vl, it diverges as L increases. 
Instead, if the gap behaves like 



A ~ exp ( — 



{Ji - Jic) a - 



(•5) 



as in the case of universality class of KT transitions (er = 
i), the Roomany-Wyld approximation for the Callen- 
Symanzik /3-function [E6|, which is defined as 



{■h) = 



1 + log I 



A L (J 4 ) 



- A L( J 4)A^ +2 (J 4 ) - 
.Ax,(J 4 )Ai + 2(J4)- 



should have the form 



Pl,L+2{Ja) ~ {Ji — Jit 



L.L+2) 



(6) 



(7) 



Fitting the form (0) to the calculated function (||) for 
each L, a is estimated as follows: <Tio,i2 = 0.38 ± 0.10, 
a-12,14 = 0.43 ±0.10 and o- 14 ,i 6 = 0.49 ± 0.10. The results 

Thus we conclude the criti- 
is characterized by the 



are consistent with a = \ 
cal behavior near J^ c for m 



universality class of the KT transition. 

Furthermore we estimate the central charge c of CFT, 
using the asymptotic form of the GS energy per site 



E(L, M) - e(m) 



7T 1 
— CD, — ^ 

6 L 2 



{L -> 00), 



(8) 



where v s is the sound velocity which is the gradient of the 
dispersion curve at the origin. The result shown in Fig. 
H suggests c = 1 with only a few percent errors for m = h 
and < Ji < 0.1. It also supports the KT transition. 

The critical exponent 77, associated with the spin cor- 
relation function in the leg direction like {SqS~) ~ 
(— l) r r~ v , can be estimated by the form of the gap 
irv s r]/L {L — > 00) 



A 



The estimated rj, shown 



in Fig. |5|, seems close to ^ around the critical point Ji c , 
rather than \ which is expected for the KT transition. 
We think there is a possible jump from r) = \ to r\ = \ 
at J^ c because the elementary excitations is expected to 
behave like the free Fermion systems (77 = |) at the edge 
of the plateau for J4 > J 4c || . The present small cluster 
analysis could not detect such a discontinuity. Another 
exponent rf defined as {SqS*) ~ cos(2fc^ r)r~ n can also 
be estimated from the L-dependence of the soft mode 



3 



gap with the momentum 2Uf = 2irm pljj . We checked 
the validity of the relation rjif = 1 around J^ c which is 
consistent with the Luttinger liquid theory ]27|] leading 
to r\ — \ in the free Fermion case. 
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FIG. 5. Estimated central charge c and exponent r\ around 
Jic. The result indicates c = 1 which is consistent with the 
KT transition, rj is close to \ rather than j. 

The spin gap at m = has already been ob- 
served in several real ladder compounds, for exam- 
ple Cu 2 (C 2 H 12 N2) 2 Cl 4 and La 6 Ca 8 Cu 24 04i @. 
The magnetization plateau, however, has not been de- 
tected at any finite magnetization. We hope some new 
ladder materials with the field-induced spin gap will be 
discovered in the near future. 

In summary the finite cluster calculation and size scal- 
ing study showed that the S = h antiferromagnetic 
spin ladder with the four-spin cyclic exchange interac- 
tion at every plaquette has the magnetization plateau at 
to = 1/2 for J4 > Ji c = 0.05 ± 0.04 and the phase tran- 
sition with respect to J4 belongs to the KT universality 
class. 
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